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SUMMARY 


This  report  considers  the  problem  of  predicting  target  coord 
nates,  given  a  temporal  and  spatial  distribution  of  shots  from  a 
multiplicity  of  targets.  Our  method  is  based  on  analyzing  multi¬ 
fold  polarity-coincidence  correlations  of  acoustic  signals  re¬ 
ceived  at  a  configuration  of  azimuth- letermining  arrays. 

We  review  the  process  of  generating  multifold  correlation 
functions,  with  emphasis  on  those  properties  which  are  of' concern 
in  azimuth  determination.  We  describe  an  azimuth  determination 
algorithm  designed  to  make  most  effective  use  of  the  information 
contained  in  the  correlation  function. 

Our  discussion  of  generalized  triangulation  points  out  the 
weaknesses  of  traditional  methods,  and  outlines  a  probabilistic 
approach  to  triangulation  that  would  yield  an  optimal  prediction 
for  a  given  set  of  azimuths.  because  this  probabilistic  approach 
would  be  extremely  demanding  in  computational  requirements,  we 
present  an  alternative  method  based  on  similar  concepts  but  yield¬ 
ing  a  closed-form  solution.  This  new  method  of  generalized  tri¬ 
angulation  is  tested  with  synthetic  data,  and  found  consistently 
more  accurate  than  the  traditional  method. 
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FOREWORD 


Thin  report  document r  progress  in  research  and  development 
of  a  signal  processing  technique  for  sound  ranging.  Some  of  the 
technical  discussion  requires  a  general  knowledge  of  previous  re¬ 
ports;  but  the  reader  unfamiliar  with  those  reports  should  still 
find  the  problems  discussed  here  both  coherent  and  generally  in¬ 
teresting. 

Overall,  this  investigation  aims  at  designing  a  data-process - 
ing  system  that  will  perform  most  accurately  and  rapidly  the  sig¬ 
nal  analysis  required  by  the  BBN  sound-ranging  system.  In  the 
previously  reported  work,  *  described  a  system  capable  of  yield¬ 
ing  correlograms  automatically  In  virtually  real  time.  These 
correlograms  represent,  .  a  1'unction  of  the  azimuthal  angle,  the 
values  of  a  fourfold  correlation  of  acoustic  signals  received  at 
the  microphones  of  an  array.  The  laboratory  sound- ranging  system 
generates  a  correlogram  approximately  every  half  second,  with  a 
resolution  of  one-half  degr* of  azimuthal  angle.  n  the  absence 
of  noise,  maxima  in  the  correlograms  should  occur  at  angles  which 
correspond,  within  ±  1/4°,  to  the  azimuths  from  which  sounds  are 
received.  In  the  current  report,  we  begin  by  assuming  many  cor¬ 
relation  functions.  Our  problem  is  to  obtain  valid  target  coordi¬ 
nates  from  such  a  sequence  of  correlograms. 
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TECHNICAL  DISCUSSION 


1.0  INTRODUCTION 

The  process  of  sound  ranging  locates  a  target  by  means  of  a 
multi-array  configuration  of  azimuth-determining  arrays.  In  pre¬ 
vious  reports , 1 * 2  *  3  we  described  how  each  array  yields,  as  a  func¬ 
tion  of  azimuthal  angle,  a  multifold  polarity-coincidence  correla¬ 
tion  (PCC)  count  of  the  signals  detected  at  its  constituent  re¬ 
ceivers.  This  report  examines  the  problem  of  determining  weighted 
sets  of  target  coordinates  from  the  information  contained  in  these 
correlation  functions.  The  procedure  to  be  followed  consists  of 
a  series  of  relatively  independent  steps.  Illustrated  in  Pig.  1 
and  outlined  below. 

Assume  that  there  exists  in  the  vicinity  of  the  sound-ranging 
system  a  number  of  targets,  each  of  which  produces  a  series  of 
shots  at  unknown  times.  .hen  the  signal  from  a  shot  impinges  on 
an  array,  it  generates  a  peak  in  the  correlation  count  obtained  at 
that  array.  As  a  first  step  then,  we  must  determine  the  direction 
of  each  shot  relative  to  each  array.  Although  one  can  approximate 
these  directions  simply  by  observing  the  angular  position  of  each 
correlation  maximum,  the  approximation  will  be  distorted  by  am¬ 
bient  noise  and  the  microphone  spectral  characteristics,  which 
give  rise  to  spurious  correlation  maxima,  and  shift  and  distort  the 
true  maxima.  In  addition,  closely  spaced  peaks  may  be  attributed 
to  either  the  distortion  of  a  single  peak  or  to  the  arrival  of  a 
number  of  shot  signals  from  almost  the  same  direction  at  almost 
the  same  time.  The  azimuth  determination  algorithm  should  take 
these  difficulties  into  account,  reject  all  spurious  azimuths,  and 
arrive  at  the  best  approximations  to  the  true  azimuths.  The  result 
of  the  azimuth  determination  algorithm  is  a  list,  for  each  array, 
of  the  azimuths  and  associated  reception  times  of  the  shot  signals. 
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The  second  step  in  the  procedure  is  to  sort  the  azimuth  lists 
and  obtain,  for  each  shot,  a  set  of  N  azimuths,  one  for  each  of 
the  N  arrays  in  the  configuration.  If  the  shots  are  amply  spaced 
in  time  and  azimuth,  this  is  a  simple  task.  But,  during  rapid-  • 
fire  sequences,  shot  signals  are  received  at  short  intervals  from 
the  same  general  direction.  Further,  the  individual  azimuth  de¬ 
terminations  are  inexact  owing  to  the  effects  of  noise  and  also 
to  the  system's  intrinsic  quantization  of  azimuthal  angle  (0.5°) 
and  reception  times  (0.5  sec).  In  some  cases,  it  may  actually  be 
impossible  to  establish  unambiguously  a  unique  set  or  sorting  of 
azimuths  for  each  detected  shot.  Through  the  use  of  heuristic 
programming,  the  system  rapidly  decides  on  a  best-possible  sort¬ 
ing  of  azimuths  in  cases  where  there  is  no  unique  correct  sorting. 
The  inherent  uncertainties  can  be  incorporated  into  a  weighting 
factor  for  each  set  of  shot  azimuths. 

In  the  ideal  case,  each  set  of  shot  azimuths  would  intersect 
at  a  single  point,  thus  specifying  uniquely  the  origin  of  each 
shot.  In  reality,  the  shot  azimuths  intersect  in  pairs  at  dif¬ 
ferent  locations;  the  derivation  of  a  pair  of  shot  coordinates 
requires  a  generalized  triangulation  algorithm.  The  specification 
of  the  shot  coordinates  by  means  of  some  average  and  variance  of 
the  shot  azimuth  intersection  points  (the  usual  approach  to  the 
problem)  has  little  basis  in  theory  and  may  lead  to  serious  dif¬ 
ficulties;  for  example,  small  errors  in  azimuth  determination  may 
induce  extremely  large  (even  infinite)  errors  in  shot  location. 

By  considering  the  bearing  lines  themselves,  rather  than  their  in¬ 
tersection  points,  one  can  derive  generalized  triangulation 
methods  -  with  a  sounder  basis  in  theory  -  that  reduce  these  dif¬ 
ficulties  . 

In  the  event  that  more  than  one  shot  emanates  from  each  tar¬ 
get,  as  a  final  step  we  must  deduce  target  coordinates  from  the 
multitude  of  shot  coordinates.  For  this  purpose,  the  target  sort¬ 
ing  algorithm  decides  how  far  apart  two  shot  locations  must  be  in 
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order  to  be  assigned  to  separate  targets,  .e  expect  that  the 
weights  generated  in  the  sorting  procedure  and  the  variances  com¬ 
puted  in  the  generalized  triangulation  procedure  will  play  a  prin¬ 
cipal  role  here,  both  in  distinguishing  the  targets  and  in  assign¬ 
ing  final  weights  to  thorn. 

In  the  entire  sound-ranging  process,  ne  must  consider  the 
influence  of  the  array  configuration  itself  on  the  accuracy  of  the 
final  target  locations.  For  a  given  error  distribution  function 
of  the  individual  azimuth  determinations,  the  target  location  er¬ 
ror  (as  a  function  of  range  and  azimuth)  will  depend  critically 
on  the  number*  of  arrays  and  on  the  manner  in  which  they  are  de¬ 
ployed.  If  we  can  obtain  an  analytical  representation  of  the  tar¬ 
get  location  error  as  a  function  of  the  multi -array  geometry,  we 
can  then  study  the  problem  of  optimizing  the  array  deployment 
using  logistic  considerations  as  constraints. 

In  the  absence  of  actual  field  data,  we  have  resorted  to  the 
use  of  synthetic  data  for  initial  study  of  tire  effectiveness  of 
the  various  algorithms.  By  assuming  shot  locations  and  times, 
meteorological  conditions,  '*rray  deployment,  ricropnone  character¬ 
istics,  and  source  and  noise  .spectra,  ;e  can  simulate  reasonably 
accurately  the  correlation  functions  which  would  be  obtained  from 
each  array  in  tire  field.  1  To  test  only  the  triangulation  scheme 
itself  we  simulated  shot  azimuths  directly  by  postulating  a  par¬ 
ticular  error  distribution  for  the  azimuth  determination  algorithm. 
Although  each  step  in  tiie  target  location  procedure  has  a  sound 
analytical  basis,  we  may  expect  that  the  use  of  such  synthetic 
data  will  be  of  great  help  in  arriving  at  the  final  form  of  the 
program. 

In  this  report,  we  limit  our  discussion  to  the  subjects  of 
azimuth  determination  and  generalized  tri angulation .  e  have  com¬ 

pleted  the  development  of  a  satisfactory  heuristic  azimuth  sorting 
algorithm,  and  preliminary  hand  analysis  with  synthetic  data  indi¬ 
cates  that  the  algorithm  yields  excellent  results  in  sorting  shots 


recorded  within  small  spreads  of  time  and  azimuth.  However,  we 
expect  that  the  algorithm  will  be  subjected  to  a  number  of  refine¬ 
ments  and  improvements  during  the  course  of  the  extensive  tests  of 
the  computer  programs  now  underway.  We  shall  therefore  postpone 
our  presentation  of  the  azimuth  sorting  algorithm  to  a  future  re¬ 
port  . 
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2.0 


AZIMUTH  DETERMINATION 


Before  proceeding  to  a  description  of  our  azimuth  determina¬ 
tion  algorithm,  we  review  briefly  the  process  by  which  the  corre¬ 
lation  functions  are  generated  and  the  difficulties  which  must  be 
overcome  in  order  to  extract  azimuths  from  them.  More  detailed 
accounts  of  these  matters  may  be  found  in  our  previous  reports. 


2.1  The  Correlation  Function:  Its  Generation  and  Characteristics 

Each  azimuth-determining  array  is  composed  of  four  T-23  mi¬ 
crophones  situated  at  the  vertices  and  center  of  an  equilateral 
triangle  with  sides  of  80  m.  (Fig.  2).  The  recorded  sound  levels 
are  analyzed  by  the  computer  in  one-second  tlme-iytes,  each  byte 
overlapping  the  preceding  byte  by  . d 6  seconds.  For  shot  signal 
durations  less  than  a  few  tenths  of  a  second,  this  procedure  en¬ 
sures  a  high  probability  of  capturing  the  entirety  of  each  signal 
in  one  or  another  byte.  Each  of  the  four  one-second  signal  sam¬ 
ples  is  converted  into  an  ordered  time-sequence  of  1000  digits: 
each  digit  is  either  +1,  0,  or  -1  depending  on  whether  the  ampli¬ 
tude  of  that  one  millisecond  sample  has  an  average  value  greater 
than  a,  between  -a  and  +a,  or  less  than  -a,  respectively.  (The 
adjustable  constant  "a"  is  the  center  clipping  level.)  The  time- 
sequences  are  analyzed  with  respect  to  azimuthal  angle  by  making 
use  of  the  azimuthal  dependence  of  the  differences  in  times  of 
arrival  of  a  signal  at  the  four  microphones.  For  each  0.5°  incre¬ 
ment  between  0°  and  180°,  we  have  derived  a  set  of  four  integers 
(incremental  time-lelay  vector)  which  specifies  the  respective 
shifts  (in  milliseconds)  of  the  four  sample  sequences  required  to 
simulate  most  closely  the  theoretically  determined  arrival  time 
differences  for  that  angle.  The  shifted  sequences  are  compared 
and  the  PCC  count  is  defined  as  the  number  of  instances  in  which 
all  four  are  +1  or  all  four  are  -1. 
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In  this  manner  v;e  generate  a  net  of  360  PCC  counts  which  may 
be  displayed  in  the  form  of  a  i  ar  graph  of  PCC  versus  azimuth. 

The  directions  of  any  shot  signals  which  reach  the  array  during 
the  one-second  tinie  byte  must  be  deduced  from  this  set  of  360  cor¬ 
relation  counts. 

The  maximum  possible  PCC  count  is  1000.  However,  for’  typi¬ 
cal  center  clipping  levels  and  shot  signal  durations,  the  maximum 
counts  will  rarely  exceed  100. 

The  major  purpose  of  the  center  clipping  is  to  reduce  tire 
deleterious  effects  of  xtraneous  noise.  In  the  absence  of  any 
shot  signal  and  with  no  center  clipping,  independent  random  noise* 
at  each  microphone  generates  an  average  PCC  count  of  125  at  each 
angle  since,  in  each  millisecond,  ther^  is  a  1/lt-  probability  that 
all  four  samples  are  positive  and  a  1/16  prolability  that  all  four 
are  negative.  ..hen  a  shot  signal  is  also  present,  this  constant 
"iCC  background  noise"  makes  the  azimuth  determination  analysis 
more  difficult.  A  center  clipping  level  "a"  -qual  to  the  rms 
noise  level  reduces  the  average  PCC  noise  contribution  by  99% . 

Of  course,  the  shot  signal  level  must  be  appreciahj_y  larger  than 
the  noise  level  in  order  that  it  not  l  e  removed  also;  .  ut  an  in¬ 
put  s ignal-to-noi se  ratio  of  one  to  two  dP  appears  to  ensure  ad‘- 
iuate  performance. 

The  followin’  Figures,  of  simulated  signal  recordings  and 
their*  associated  PCC  bar  graphs  illustrate  th*  important  aspects 
of  the  correlation  function.  In  Fig.  j,  we  reproduce  the  record¬ 
ings  whicli  would  1  e  obtained  at  the  four  microphones  of  an  array 
under  the  following  ideal  conditions:  (1)  a  typical  "H-wave"  pres¬ 
sure  pulse  incident  on  the  array  at  an  azimuthal  angle  of  65°, 

(2)  no  background  noise,  (3)  microphones  having  infinite  spectral 
response,  (A)  no  center  clipping.  The  corresponding  PCC  bar  graph 
has  a  single,  well  defined  maximum  at  the  correct  azimuth  and  is 
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SIMULATED  NOISELESS  SIGNALS  AND  PCC  FUNCTION 
PRESSURE  PULSE  AT  65° 


zero  everywhere  else.  The  maximum  PCC  count  of  only  fifteen*  in 
indicative  of  the  extremely  short  duration  of  the  shot  signal. 

Figure  ^  shows  the  effect  of  having  two  such  H-wave  pressure 
pulses  incident  at  the  same  time  but  at  different  angles  (65°  and 
7f.5°)>  all  other  conditions  the  same.  The  correlation  function 
again  has  unique  maxima  at  the  correct  azimuths. 

Figure  5  duplicates  the  conditions  in  Fig.  3>  except  that  we 
have  simulated  the  response  of  actual  T-23  microphones  to  the 
pulse.  Acting  essentially  as  narrowband  filters,  the  microphone; 
convert  the  pulse  into  a  temporally  extended  tone  burst .  The 
principal  effect  on  the  correlation  function  is  that  the  maximum 
PCC  count  is  increased  (from  13  to  ><7)  as  a  result  of  the  longer 
duration  of  the  recorded  signals. 

In  Fig.  6,  we  have  the  T-23  responses  to  the  same  two  simul¬ 
taneous  pulses  as  in  Fig.  l\.  Although  the  responses  appear  to  be 
from  a  single  shot,  the  correlogram  shows  the  ability  of  the  sys¬ 
tem  to  clearly  separate  the  two  signals. 

in  consideration  of  the  influence  of  random  noise  on  the  cor¬ 
relation  function,  Fig.  7  shows  the  effect  of  superimposing  inde- 
pendent  noise  on  the  pressure  pulses  of  Fig.  3-  In  producing  the 
correlogram,  the  center  clipping  level  "a"  was  set  so  that  the 
noise  observed  on  the  recordings,  ithout  the  signal,  exceeded  "  " 
approximately  57  of  the  time.  Comparing  the  correlograms  of 
Figs .  3  and  7,  we  note  the  following  consequences  of  the  introduc¬ 
tion  of  noise:  (1)  a  random  modulation  of  the  correlation  peak, 
which  in  turn  distorts  the  angular  location  of  the  maximum  PCC 
count,  (2)  small  localized  correlation  counts  at  angles  far  re¬ 
moved  fror.  the  signal  maximum  (the  correlation  response  for  angles 


*The  bar  graphs  in  these  Figures  represent  the  PCC  counts  normal¬ 
ized  to  the  largest  count;  the  actual  number  of  counts  corre- 
■  ponding  to  each  bar  is  indicated  by  the  number  at  its  lase. 
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SIMULATED  NOISELESS  SIC-NALS  AND  PCC  FUNCTION: 

TWO  SIMULTANEOUS  PRESSURE  PULSES  AT  65°  AND  76.5 


SIMULATED  NOISELESS  T-23  RESPONSES  AND  PCC  FUNCTION 
PRESSURE  PULSE  AT  65° 
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SIMULATED  SIGNALS, WITH  INDEPENDENT  NOISE, AND  PCC  FUNCTION 
PRESSURE  PULSE  AT  65°,  INPUT  THRESHOLD  SENSITIVITY  AT  .05 


outside  the  range  illustrated  was  no  higher  than  the  highest  noise 
response  in  the  displayed  portion)  and,  (3)  the  center  clipping 
used  to  reduce  the  noise  response  also  reduces  the  correlation  re¬ 
sponse  of  the  signal  itself  (in  this  case,  from  15  to  6). 

'.he  preceding  examples  illustrate  the  principal  difficulties 
which  must  be  overcome  by  a  satisfactory  azimuth  determination  al¬ 
gorithm;  it  must  (1)  reject  all  secondary  correlations  due  to 
noise,  (2)  recognize  and  accept  all  correlations  due  to  actual 
shot  signals,  (3)  decide  whether  overlapping  correlation  peaks 
represent  distortion  of  a  single  peak  or  the  presence  of  more  than 
one  shot  signal  and,  (h)  make  best-possible  estimates  of  the  azi¬ 
muths  of  all  shot  signals. 


2.2  The  Azimuth  Determination  Algorithm 

The  following  subsections  describe  elements  of  the  algorithm 
as  implemented  in  our  present  ystem;  we  have  also  suggested  pos¬ 
sible  improvements  for  a  future  system.  The  algorithm  involves  a 
number  of  constants  whose  values  we  have  not  specified.  The  opti¬ 
mal  choices  for  these  constants  must  await  the  exercising  of  the 
program  with  synthetic  data  or,  where  possible,  with  actual  field 
data. 


2.2.1  Correlation  threshold  and  maximum  test 

This  test  is  designed  to  eliminate  the  small  local  PCC  counts 
caused  by  noise  and  to  identify  any  present  signal  correlation 
peaks.  We  define  the  mean  correlation  M  as  the  average  of  the 
360  PCC  counts.  The  minimal  signal  correlation  S  is  defined  as 
a  linear  function  bM  +  c  of  the  mean.  If  the  largest  count  is 
less  than  S,  we  assume  that  no  signal  is  present  and  the  analysis 
is  terminated.  The  correlation  threshold  T  is  defined  as  another 
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linear  function  dM  +  e  of  the  mean,  ouch  that  only  a  small  frac¬ 
tion  (perhaps  10%)  of  all  the  nonzero  counts  are  greater  than  T. 

The  threshold  is  subtracted  from  each  count,  thus  removing  most, 
arid  possibly  all,  of  the  noise  contribution. 

2.2.2  Group  weight  test 

A  correlation  group  is  defined  as  a  set  of  three  or  more  con¬ 
secutive  nonzero  PCC  counts.  All  counts  which  do  not  belong  to 
any  group  are  eliminated.  he  weight  of  a  group  is  the  sum  of  its 
component  counts.  The  average  weight  of  a  group  is  its  weight 
divided  by  the  number  of  counts  it  contains.  Any  group  whose  aver¬ 
age  weight  is  less  than  a  predetermined  minimum  weight  V/  is  elimi¬ 
nated.  This  test  should  remove  any  vestigal  noise  contributions. 


2.2.3  Group  azimuths  and  group  overlap  test 

Starting  at  the  beginning  of  each  group,  we  add  the  counts  in 
the  group  until  the  sum  most  nearly  equals  half  the  group  weight. 
The  azimuth  of  the  last  count  in  the  sum  is  the  group  azimuth. 

Thus,  the  azimuth  associated  with  each  group  is  the  azimuth  which 
bisects  the  "area"  f  the  group  (which  is  not  necessarily  the  azi¬ 
muth  having  the  largest  PCC  count  in  the  group).  The  most  prob¬ 
able  position  of  this  half-area  azimuth  is  the  position  of  the 
original  (noiseless)  peak,  since  the  random  noise  adds,  on  the 
average,  as  many  counts  to  one  sid^  of  the  original  peak  as  to  the 
other.  Thus,  the  above  procedure  should  provide  an  effective  means 
of  allowing  for  noise-induced  distortion  of  the  correlation  peaks. 

If  either  (a)  the  PCC  count  at  the  group  azimuth  is  the  maxi¬ 
mum  count  in  the  group,  r  (b)  the  maximum  count  in  the  group  ex¬ 
ceeds  the  group-azimuth  count  by  less  than  a  predetermined  overlap 
test  constant  G,  then  the  original  group  azimuth  is  retained. 
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However,  if  the  maximum  count  in  the  group  exceeds  the  group- 
azimuth  count  by  more  than  G,  we  subtract  the  group-azimuth  count 
from  each  count  in  the  group  and  repeat  the  analysis  (i.e.,  find 
group  azimuth,  do  overlap  test)  for  each  of  the  new  groups  thus 
generated.  In  Fig.  8,  we  illustrate  the  overlap  test  with  a  typi¬ 
cal  example  of  overlapping  correlation  groups. 

We  may  expect  that  this  analysis  will  fail  if  two  signal 
groups  overlap  to  such  an  extent  that  their  "sideskirts"  sum  to 
a  larger  count  than  either  of  their  maxima.  It  .’ould  then  be  ex¬ 
tremely  difficult  to  detect  the  presence  of  two  separate  shot  sig¬ 
nals.  This  problem  might  be  circumvented  under  certain  conditions 
by  making  use  of  the  shot  azimuths  determined  at  the  other  arrays 
of  the  configuration  during  the  same  time  byte.  If  the  two  sig¬ 
nal  groups  at  the  other  arrays  have  sufficient  angular  separation 
to  be  resolved,  it  would  then  be  clear  that  there  was  a  "missing" 
azimuth  at  one  array;  the  single  azimuth  found  at  that  array  would 
be  counted  twice  in  the  azimuth  sorting  process. 


2.2.4  Use  of  overlapping  azimuth  intervals 

In  the  algorithm  described  above,  we  analyze  all  counts  in 
the  entire  180°  azimuth  span  at  once;  but  we  can  reduce  the  com¬ 
puter  storage  requirements  by  dividing  the  correlation  function  in¬ 
to  a  number  of  azimuthal  segments  and  exercising  the  algorithm  on 
each  segment  separately.  To  ensure  that  individual  groups  are  not 
split  up,  the  segments  should  overlap  by  an  amount  that  is  large 
in  comparison  to  the  maximum  expected  azimuthal  spread  of  a  signal 
group.  We  currently  use  five  50°  segments  with  17-5°  overlaps. 


PCC  COUNT  ^  ^  PCC  COUNT 


AZIM 


PERFORM  OVERLAP  TE 


SUBTRACT  GROUP  - AZ I 


1 0  )  =  0 
i  3  )  =  0 


3.0  GENERALIZED  TR I ANGU LAT ION 


When  we  have  a  set  of  shot  azimuths  generated  by  several 
shots  at  different  targets  and  no  a  priori  knowledge  of  the  cor¬ 
respondence  of  shots  and  azimuths,  we  must  first  select  subsets 

Ni 

of  shot  azimuths  corresponding  to  individual  shots.  By  azimuth 
sorting,  we  obtain  sets  of  N  shot  azimuths  (assuming  an  N-array 
configuration),  each  set  corresponding  to  a  unique  shot.  We  must 
then  consider  the  most* appropriate  generalized  triangulation  pro¬ 
cedure  for  deducing  a  shot  location  given  a  set  of  N  intersecting 
bearing  lines  . 

Traditional-  sorting  methods  tend  to  rely  to  some  degree  on 
finding  azimuths  whose  pairwise  intersection  points  are  grouped 
in  small  clusters,  .'uch  sorting  methods  involve  implicit  assump¬ 
tions  that  the  subsequent  generalized  triangulation  procedure 
should  be  based  primarily  on  the  positions  of  the  intercepts,  and 
that  the  closeness  of  intercepts  is  a  valid  criterion  for  the  re¬ 
liability  of  a  target  prediction.  In  this  Section,  we  shall  show 
that  these  assumptions  do  not  necessarily  lead  to  realistic  tar¬ 
get  locations.  Instead,  a  proper  generalized  triangulation  scheme 
should  emphasize  the  relationship  of  the  shot  prediction  to  the 
bearing  lines  themselves,  rather  than  to  their  intersection  points. 
In  a  future  report,  we  will  explore  this  concept  in  connection 
with  the  sorting  procedure. 


3.1  Generalized  Tri angulation  Using  Azimuth  Intersections 

Consider  first  the  simplest  case  of  just  two  intersecting 
bearing  lines.  The  point  of  intersection  is  then  the  only  reason¬ 
able  shot  location  prediction.  The  displacement  of  this  predic¬ 
tion  from  the  true  shot  location  depends  on  the  accuracy  of  the 
bearings  and  on  the  angle  at  which  they  intersect .  uppose  that 
the  bearings  intersect  at  a  small  angle  <t> ,  that  the  uncertainty  of 
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each  bearing  is  an  angle  A<p,  that  the  range  from  the  center  of  the 
array  configuration  to  the  intersection  point  is  R,  and  that  R  is 
much  greater  than  the  interarray  spacing.  Then  there  is  an  un¬ 
certainty  of  approximately  2RA$/$  in  the  range  to  the  shot,  and  an 
uncertainty  of  approxi  mat'  ly  RA in  a  direction  perpendicular  to 
the  direction  of  the  range.  For  shot  locations  whose  bearings  sub¬ 
tend  small  angles  with  respect  to  the  line  joining  the  two  arrays, 
the  range  uncertainty  is  a  large  fraction  of  the  range  itself. 

/.e  might  intuitively  expect  that,  if  we  obtain  a  third  bear¬ 
ing,  the  triangle  determined  by  the  addition  of  the  third  bearing 
would  have  a  high  probability  of  enclosing  the  true  shot  location. 
To  understand  why  this  is  not  in  fact  tie  case,  let  us  assume  that 
the  bearing  errors  ar-j  distributed  symmetrically;  i .  e . ,  the  true 
shot  location  is  eaually  lik<-ly  to  he  on  cither  side  of  any  bear¬ 
ing  line.  it  follows  then  that  the  prol  al  illty  of  the  shot  being 
in  any  one  of  the  four  sectors  created  o.v  the  intersection  of  two 
bearing  liras  i.  one-fourth.  third  hearing  line  could  determine 
a  triangle  by  closing  off  either  of  two  opposing  sectors.  The 
triangle  formed  by  c Losing  off  cither  sector  would  have  a  1/8 
probability  of  enclosing  the  shot.  The  total  probability  that 
the  shot  i s  within  the  triangle  is  therefore  one-fourth.  With 
only  .  1  "  confi  i^nc-  that  tic  shot  is  in  th<  triangle,  we  conclude 
that  the  triangle  -»-n<  rated  by  the  intersections  of  three  bearing 
lines  is  a  questl<  nabl e  improvement  over  th<^  single  point  gener¬ 
ated  by  t!  e  inter,  •  r*  ion  of  i  w  'nearing  liter. . 

■.'hen  an  even  larger  number  of  bearings  are  available,  the 
intersection  points  may  be  aistributed  ovt  r  a  large  area  and  will 
not  usually  a  Single  closed  polygon.  An  intuitive  approach 

might  lead  us,  to  specify  the  Shot  location  by  means  of  some  aver¬ 
age  of  all  the  intersection  points.  in  the  next  subsection,  we 
inscribe  an  alg  ri  *  t  ased  on  tills  ap:  roach  and  show  some  of 
the  difficultie  a.  sociated  with  it. 
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3.1.1  Method  of  the  mean  of  the  intersections  (MI) 

A  set  of  N  shot  azimuths  in  general  defines  N(N-l)/2  inter¬ 
section  points;  the  number  of  intersections  is  smaller  if  some  of 
the  azimuths  are  parallel.  Given  the  position  of  each  azimuth- 
determining  array,  the  x  and  y  coordinates  of  the  intersection 
points  relative  to  some  fixed  origin  can  be  calculated  by  means 
of  standard  trigonometric  formulas.  The  shot  location  coordinates 
are  obtained  by  averaging  separately  the  x  and  y  coordinates  of 
all  the  intersection  points.  We  henceforth  refer  to  this  algorithm 
as  the  MI  method. 

3.1.2  Inadequacies  of  the  MI  method 

The  most  obvious  fault  of  the  MI  method  is  that  it  fails  com¬ 
pletely  if  any  azimuths  are  parallel:  the  predicted  shot  location 
is  then  at  infinity. 

Another  drawback  is  that  the  prediction  may  be  very  sensitive 
to  small  errors  in  the  individual  azimuth  determinations.  Con¬ 
sider,  for  example,  the  intersections  of  the  four  bearing  lines 
shown  in  Fig.  9.  The  solid  lines  create  six  intersections  whose 
mean  is  shown  ty  the  dot  enclosed  by  a  solid  circle.  The  dashed 
lines  represent  an  alternative  set  of  bearing  lines  chosen  at  ran¬ 
dom  under  the  condition  that  their  directions  differ  from  the 
original  directions  1  y  no  more  than  a  small  assumed  angular  uncer¬ 
tainty.  The  alternate  tearing  lines  form  intersections  that  shift 
the  prediction  to  the  point  indicated  by  the  dot  enclosed  by  a 
dashed  circle.  Assume  that  the  original  prediction  is  close  to 
the  true  shot  location.  Although  two. of  the  alternate  bearing 
lines  pass  closer  to  the  shot  than  the  corresponding  solid  lines, 
the  alternate  prediction  is  nonetheless  poorer  than  the  original 
predict  ion . 
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A  common  method  of  reducing;  these  difficulties  is  simply  to 
discard  (or  assign  small  weight  to)  intersections  which  are  far 
removed  from  the  mean.  his  is  a  rather  artificial  procedure  and, 
moreover,  it  offers  no  guarantee  that  the  prediction  obtained  from 
the  surviving;  azimuths  will  be  better  than  the  original  prediction. 
Furthermore,  if  w<-  have  only  three  bearing  lines,  there  is  no  ob¬ 
vious  criterion  for  determining^^iri-ch  intersection  is  the  "bad" 
one . 


3.2  The  Probabilistic  Approach 

The  preceding  discussion  leads  to  the  conclusion  that  the  Ml 
triangulation  method  is  not  based  on  sound  analytical  concepts. 
Indeed,  thus  far  little  fundamental  research  lias,  teen  devoted  to 
the  problem  of  determining  an  optimal  generalized  trihngulation 
algorithm.  We  outline  below  an  extremely  general  approach  to  this 
problem,  the  implications  of  which  we  hope  to  develop  more  fully 
in  our  future  work. 

Our  main  point  is  that  the  prediction  of  a  shot  location  is 
essentially  a  statistical  process.  ..e  are  presented  with  a  set 
of  bearing  lin^s  and  associated  statistical  uncertainties;  our 
goal  should  be  to  determine  that  point  on  the  plane  which  has  the 
greatest  probability  of  coinciding  with  the  shot  location.  To 
accomplish  this,  we  may  transform  the  probalility  functions  of 
the  individual  azimuth  uncertainties  into  a  probability  function 
for  points  on  the  plane.  The  resulting  two-dimensional  probabil¬ 
ity  function  would  contain  a  large  amount  of  useful  information. 

Its  maximum  would,  of  course,  determine  the  best  shot  location 
prediction.  The  height  at  the  maximum  point,  the  steepness  of  the 
probability  surface  in  the  vicinity  of  the  maximum  point,  and  the 
mean  level  of  the  probability  function  outride  the  region  of  the 
maximum  would  be  some  measure  of  the  reliability  of  the  prediction. 
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A  simple  way  to  visualise  the  process  is  to  imagine  a  set  of 
translucent  overlay®,  on  each  of  which  is  imprinted  a  variable  ' 
density  optical  image  of  a  bearing  line  and  its  associated  proi - 
ability  density.  The  image  is  most  dense  along  the  bearing  line. 
The  density  decreases  for  azimuths  on  either  side  of  the  bearing 
line  and  approaches  zer  (i.  . ,  transparency)  for  azimuths  further 
from  the  bearing  line  than  the  assumed  uncertainty  of  the  azimuth 
determination.  The  overlays  are  placed  on  top  of  one  another  on 
a  target  map  in  such  a  way  that  each  bearing  line  passes  through 
the  appropriate  array  posi*  ion  at  the  correct  angle.  The  point 
on  the  map  at  which  the  combined  overlays  are  most  opaque  would 
then  be  the  most  prot  a:  lr  shot  location. 

The  analytical  Implementation  of  the  process  requires  an  as¬ 
sumption  as  to  the  statistics  of  the  azimuth  uncertainties.  e 
could  assume  a  Gaussian  distribution  of  the  deviation  of  an  azi¬ 
muth  from  the  rm*an,  ..  ith  a  variance  proportional  to  the  expected 
uncertainty.  We  must  then  find  the  geometrical  transformation 
between  x,y  coordinates  r<  latlve  to  some  fixed  origin,  and  range 
and  azimuthal  deviation  relative  to  the  array.  The  Jacobian  of 
this,  transformation  is  required  to  perform  the  transformation  of 
probability  densities  to  x,y  coordinates.  1  "n  this  is  done  for 
each  array  in  the  configuration,  tin1  final  two-dimensional  prob¬ 
ability  function  is  calculate  1  as  the  product  of  the  individual 
functions. . 

because  this  method  is  based  on  the  statistics  of  the  azimuth 
directions  themselves,  and  not  on  the  statistics  of  the  azimuth 
intersection  points,  .-.e  would  expect  that  the  prediction  obtained 
would  be  quite  insensitive  to  small  Individual  azimuth  uncertain¬ 
ties. 

The  principal  drawback  of  the  general  method  is  that  although 
quite  simple  in  concept,  it  would  :  <_•  extremely  demanding  in  compu¬ 
tational  requirements,  and  thus ,  beyond  the  scope  of  the  present 


system.  In  the  next  Section,  we  outline  an  alternative  method 
which  also  relies  on  the  bearing  lines  themselves  but  is  more 
feasible  in  practice. 


3.3  Generalized  T ri an g u 1  at i on  Using  Azimuth  Displacements 

In  lieu  of  the  rigorous  probabilistic  approach  described 
above,  we  seek  a  generalized  triangulation  procedure  which  satis¬ 
fies  the  following  criteria: 

1.  .tie  procedure  should  be  based  on  some  relation  involving  the 
bearing  lines  themselves,  rather  than  their  intersection 
points  . 

.  The  shot  prediction  should  be  relatively  insensitive  to  small 
deviations  of  the  bearing  line  directions. 

* 

3.  The  method  should  be  reasonable  in  computational  requirements; 
preferably,  the  coordinates  of  the  shot  prediction  should  lie 
expressible  as  a  closed  analytical  formula,  rather  than  as  the 
result  of  a  search  procedure. 

A  method  which  seems  to  satisfy  all  of  these  criteria  is  based  on 
the  perpendicular  distances  from  an  assumed  shot  location  to  all 
of  the  bearing  lines. 


3.3.1  Method  of  the  minimum  mean-square  displacement  (MD) 


The  MD  prediction  in'  defined  as  that  point  which  minimizes 
the  mean-square  perpendicular  distance  to  all  bearing  lines.  In 
Appendix  A,  we  derive  simple  closed  formulas  expressing  the  coor¬ 
dinates  of  the  MD  prediction  in  terms  of  the  coordinates  of  the 
N  arrays  and  the  azimuths  of  the  N  bearing  lines. 


The  actual  value  of  the 
prediction  from  the  bearings 


root-mean- -quare  displacement  of  the 
might  be  a  suitable  measure  of  the 


25 


uncertainty  to  be  assigned  to  the  prediction.  However,  other  mea¬ 
sures  of  uncertainty  should  also  be  considered.  A  valid  uncer¬ 
tainty  measure  should  be  well  correlated  with  the  actual  distance 
from  the  prediction  to  the  true  shot  location. 

n  the  following  subsection,  we  compare  the  relative  sensi¬ 
tivities  of  the  MI  and  MD  methods  to  small  deviations  of  the  bear¬ 
ing  line  directions. 

3.4  Qualitative  Comparison  of  MD  and  MI  Predictions 

We  may  illustrate  the  major  ivantages  of  the  MD  method  by 
means  of  some  simple  qualitative  examples.  First,  we  consider  a 
linear  deployment  of  three  arrays  with  an  interarray  separation 
D,  and  a  shot  originating  on  the  central  axis  of  the  configuration 
at  a  distance  HD  (Fig.  10).  ith  a  perfect  azimuth  determination 
routine,  the  three  bearing  lines  would  intersect  at  the  shot  loca¬ 
tion  and  either  the  MI  or  MD  method  would  yield  a  perfect  predic¬ 
tion  (Fig.  10a).  We  now  assume  progressively  larger  errors  in 
the  bearing  line  associated  with  the  leftmost  array  so  that  it 
gradually  becomes  parallel  to  the  center  bearing  line.  Figure  10b 
shows  that  if  the  perturbation  is  small,  the  predictions  of  the 
two  methods  almost  coincide.  As  the  enclosed.,  triangle  becomes 
more  elongated,  as  in  Fig.  10c,  the  MI  prediction  moves  out  from 
the  true  shot  location  much  more  rapidly  than  does  the  MD  predic¬ 
tion.  when  the  bearing  lines  become  parallel,  as  in  Fig.  lOd,  the 
MI  prediction  moves  out  to  infinity,  while  the  MD  method  still 
produces  a  finite  prediction  (although  at  considerable  distance 
from  the  shot,  displaced  by  approximately  a  50%  increase  in  range). 

We  can  construct  another1  example  for  which  the  MI  method 
fails  completely  while  the  Ml  method  actually  yields  a  perfect 
prediction.  We  consider1  a  linear  four-nrray  configuration  with 
the  shot  again  on  the  central  axis  (Fig.  11)  and  study  the  effect 
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of  errors  which  tend  to  make  the  two  inner  bearing  lines  parallel. 
The  MI  prediction  again  rapidly  diverges  as  the  bearing  lines  ap¬ 
proach  parallelism.  The  MD  prediction  moves  only  a  small  maximum 
distance  from  the  shot  location  and  then  moves  back  to  coincidence 
with  the  shot  location  in  the  limiting  case.  The  example  admit¬ 
tedly  has  an  extremely  unlikely  symmetry;  nevertheless,  it  dra¬ 
matically  illustrates  the  general  tendency  of  the  MU  method  to 
retain  a  reasonable  prediction  even  when  the  MI  method  fails. 

These  examples  illustrate  the  following  important  general 
result:  an  MD  prediction  exists  as  long  as  at  least  two  bearing 
lines  are  not  parallel.  This  conclusion  also  follows  directly 
from  the  analytical  formulas  in  Appendix  A.  In  contrast,  the  MI 
method  yields  a  prediction  at  infinity  with  the  presence  of  any 
parallel  bearing  lines. 

Another  property  of  the  MD  method  is  that  two  parallel  bear¬ 
ings  may  be  replaced  by  a  single  bearing  originating  at  a  virtual 
array  located  midway  between  their  respective  arrays  —  and  the 
prediction  will  remain  unaltered.  For" example ,  in  Fig.  lOd  the 
prediction  is  at  the  intersection  of  the  bearing  from  the  right¬ 
most  array  and  a  bearing,  parallel  to  the  two  parallel  bearings, 
emanating  from  a  virtual  array  midway  between  the  leftmost  array 
and  the  center  array.  Figure  lid  can  be  described  in  a  similar 
manner. 


3.5  Quantitative  Comparison  of  MD  and  MI  Predictions 

The  generation  of  the  synthetic  data  used  to  test  our  gen¬ 
eralized  triangulation  algorithm  is  described  in  Appendix  D.  /.e 
assumed  a  four-array  configuration  consisting  of  arrays  at  the 
vertices  and  center  of  an  equilateral  triangle  of  side  875  m. 
(Fig.  -1).  We  postulated  a  set  of  19  targets  at  various  ranges 
and  azimuths  from  the  center  of  the  configuration  (Fig.  E-2). 
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The  ranges  vary  between  1.5  and  ^.0  km.  We  further  postulated  that 
a  total  of  37  shots  originate  at  the  targets  at  specified  times 
(Table  B-I).  We  calculated  the  precise  azimuth  of'  each  shot  rela¬ 
tive  to  the  center  of  each  array  (Table  B-  I)  and  then  polluted  the 
data  to  simulate  the  azimuthal  errors  to  be  expected  in  a  real  sys¬ 
tem  (Table  B-IIJ).  .e  also  computed  the  arrival  times  (rounded  off 
to  0.5  sec)  of  each  shot  at  each  array. 

In  this  subsection,  we  study  the  results  of  applying  both  the 
MI  and  MD  methods  to  each  of  the  37  sets  of  shot  azimuths  in 
Table  B-TII.  (V/e  realize  that  in  doing  this,  we  are  temporarily 
ignoring  the  azimuth  sorting  problem;  we  know  that  each  of  the  37 
shot  azimuths  sets  correponds  to  a  unique  shot  because  of  the  way 
in  which  they  were  gen'  rated.) 

The  results  of  the  calculations  are  presented  in  Table  I. 

For  each  shot,  we  give  the  difference  in  range  and  azimuth  (rela¬ 
tive  to  the  centei*  of  the  configuration)  between  the  actual  and 
predicted  shot  locations.  We  also  give  the  radial  error,  .hich  is 
the  total  distance  of  the  prediction  from  the  true  shot  location. 

In  Fig.  1?,  wo  plot  scale  map.  of  the  vicinity  of  target  1,  show¬ 
ing  the  converging  shot  azimuths  and  the  locations  of  the  MI  and 
MD  predictions  for  Shots  4,  8,  and  30.  Target  1  is  at  a  range  of 
5000  in.,  so  that  the  position  of  the  center  of  the  array  config¬ 
uration  would  he  11  in.  to  the  left  of  the  target  on  these  maps. 

The  arx’ay  configuration  would  be  about  2  in.  on  a  side. 

The  >'vror  quantity  of  greatest  Interest  in  sound  ranging  is 
the  normalized  radial  error*  defined  as  the  ratio  of  the  radial 
error  to  the  actual  range  of  the  shot.  From  the  data  of  Table  I, 
we  calculated  an  av^rag^  normalized  radial  error  of  3-9%  for  the 
37  MD  predictions.  In  considering  the  MI  predictions,  we  are 
faced  with  the  difficulty  that  two  predictions  are  infinite  due  to 
the  existence  of  parallel  bearing  lines.  Ignoring  these  tiro  fail¬ 
ures  and  averaging  the  other  35  MI  predictions ,  v.e  obtained  an 
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TABLE  I.  Comparison  of  MD  and  MI  shot  location  predictions. 


Range 

(km) 

Radial  Error 

Range  Error 

Azimuthal  Error 

■■ 

MD 

(m) 

MI 

(m) 

MD 

(dec) 

MI 

(dec) 

1.5 

6 

131 

-  7 

131 

0.0 

.  2 

1.5 

18 

15 

17 

-  6 

.2 

.5 

1.5 

14 

32 

1 

-  30 

-  -5 

- 

.4 

5.0 

311 

462 

310 

462 

•  3 

.  1 

4.  0 

152 

1113 

-152 

-1110 

0.0 

-1 

.3 

4.0 

235 

225 

-235 

-  223 

0.0 

- 

.4 

4.0 

122 

545 

-121 

-  545 

0.0 

.  1 

5.0 

390 

608 

-390 

-  608 

.3 

.1 

4.o 

160 

572 

-160 

-  572 

0.0 

- 

.4 

3.0 

14  5 

INFINITE 

-145 

INFINITE 

0.0 

.5 

3.0 

154 

305 

-154 

-  305 

0.0 

- 

.1 

3.0 

46 

183 

-  42 

-  182 

_  #  'J 

- 

•  3 

3.!' 

64 

109 

63 

-  108 

0.0 

.3 

3.0 

14  4 

74 

-14  4 

-  74 

.3 

- 

.2 

3.0 

155 

86 

154 

86 

0.0 

- 

.  1 

3.3 

219 

339 

-210 

-  339 

0.0 

.1 

3.4 

3 

37 

2 

-  36 

.  1 

.2 

3-5 

311 

558 

311 

558 

.  1 

- 

.2 

3.6 

75 

120 

73 

-  120 

-  .3 

- 

.0 

3.7 

207 

1409 

204 

1408 

• 

- 

.6 

5.0 

142 

107 

136 

-  99 

J 

•  J 

.5 

4  .0 

160 

572 

-160 

-  572 

0.0 

- 

.4 

4.0 

78 

42  4 

-  7° 

424 

0. 

- 

.2 

4.0 

122 

545 

-121 

-  545 

0.0 

.1 

2.0 

14 

164 

13 

164 

-  .3 

- 

.1 

4. 

152 

555 

-153 

-  555 

0.0 

- 

5.0 

506 

48 

505 

-  24 

-  .2 

- 

•  5 

4.0 

76 

39 

-67 

5 

.5 

.6 

4.0 

76 

INFINITE 

-  75 

INFINITE 

-  .5 

3 

•  5 

5.0 

88 

245 

84 

24  5 

•  J 

.  1 

5.0 

224 

235 

222 

-  2  35 

.3 

- 

.0 

4.0 

152 

1113 

-152 

-  1110 

0.0 

-1 

.3 

4.0 

334 

216 

334 

-  219 

0.0 

- 

.2 

4.0 

186 

273 

-186 

273 

0.0 

- 

.2 

2.0 

61 

112 

-  61 

112 

.2 

.  0 

4.0 

262 

66 

261 

64 

0.0 

.2 

2.0 

15 

28 

-  1 

-  23 

-  •  3 

- 

:5 

average  normalized  radial  error  of  8.9?  for  the  MI  method.  In 
summary,  we  find  an  average  radial  error  of  about  ^0  m.  \  er  km. 
for  the  MD  method  and  more  than  90  m.  pei*  km.  for  the  MI  method. 

The  MD  [-rediction  is  superior  to  the  MI  prediction  in  fully 
28  of  the  37  cases. 

..e  are  currently  carrying  out  a  similar  analysis  using  syn¬ 
thetic  data  generated  with  an  assumed  set  of  targets  having  ranges 
of  up  to  several  dozen  kilometers.  We  expect  that  the  results  will 
not  differ  markedly  from  those  obtained  here. 
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4.0  RECOMMENDATIONS 


We  have  outlined  the  general  procedure  for  predicting  target 
locations  by  means  of  multifold  correlations  of  received  acoustic 
signals.  We  have  considered  in  detail  the  steps  in  the  procedure 
concerned  with  azimuth  Jetermination  and  generalized  triangula¬ 
tion.  The  Sorting  problem  will  lo  considered  in  a  future  report. 

In  regard  to  azimuth  determination,  we  recommend  that  further 
refinements  of  our  algorithm  be  examined  when  the  result:;  of  field 
tests  are  available.  In  particular,  such  field  data  should  enable 
us  to  determine  the  optimal  values  of  the  variou.  test  constants 
involved  in  the  algorithm. 

In  regard  to  generalized  trl angulation ,  we  recommend  that  the 
method  of  minimum  mean-  quar^  displacement  (MI  )  I  used  in  future 
studies.  In  addition,  the  mor<>  rigorous  probabilistic  approach 
should  be  studied  further  in  or  ler  to  doterm'in'  whether  it  might 
be  Implemented  In  a  tract  leal  system. 
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APPENDIX  A 


Derivation  of  Formulas  for  the  MD  Method 
of  Generalized  T ri an gul at i on 

We  wish  to  find  expressions  for  the  coordinates  of  that  point 
v/hose  mean-  •quare  disf  lacement  from  a  set  of  N  bearing  lines  is 
minimal.  Figure  A-l  illustrates  the  geometry  of  the  arrays,  bear¬ 
ings,  target,  and  coordinate  origin.  The  directions  of  the  x  and 
y  axes  are  specified  by  the  unit  vectors  i  and  j,  respectively. 

The  locations  of  the  arrays  relative  to  the  origin  are  specified 
by  the  vectors  (ot  =  l,  ,...,11).  The  direction  of  the  bearing  line 
determined  at  the  ath  array  is  specified  by  the  unit  vector  n  , 
which  makes  an  angle  <p  with  the  y  axis. 

The  point  to  b  ■  determined  is  specified  by  the  variable  vec¬ 
tor  R,  relative  to  t  he  origin,  and  by  the  variable  vector  xa>  rela¬ 
tive  to  the  ath  array.  'Ine  displacement  of  the  point  from  the  ath 

bearing  line  is  given  by  the  magnitude  of  the  vector  t  .  The  dis- 

0 1 

tance  from  the  ath  array  to  the  foot  of  the  vector  Qa  is  given  by 

d  . 

a 

fince  ria  is  a  unit  vector  perpendicular  to  t.  ,  we  have 


n 

a 


1 


a 


e  =  0  . 

-  a 


(A-l) 


Since  v  can  be  expressed  either  as  R  -r  or  as  e  +  afi  ,  we  ob- 
~a  ~~a  -  a  a  a 

tain 


e 

-  a 


R 


r 

~a 


A 

d  n 


a  a 


( A-Z  ) 


Taking  th«  scalar  product  of  Eq.  A-2  with  the  unit  vector  ,  and 
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ct:  o 


applying  Kq .  A-l,  we  obtain  an  expression  for  dQ: 


da  =  •  (A~3) 

If  we  substitute  this  expression  back  into  Eq.  A-?,  and  then  form 
the  scalar  product  of  the  resulting  equation  with  itself,  we  ar¬ 
rive  at  an  expression  for  the  squared  magnitude  of  the  displace¬ 
ment  vector: 


[n 


a 


(0-ra)l2 


(A-U) 


The  mean-square  displacement  of  the  assumed  shot  location  from 
the  N  bearing  lines  is  thus 


<  |  H-r*  I2  >-<  [n  •  ( R-r )  1  2  > 


(A-5) 


(The  angular  brackets  denote  an  average  over  the  N  arrays  of  the 
expression  they  enclose;  for  convenience,  we  delete  all  subscripts 
within  the  brackets.) 

V/e  now  assume  that  the  function  o2  describes  a  concave  sur¬ 
face  as  a  function  of  the  variable  position  of  the  vector  R.  To 
find  the  minimum  of  this  surface,  we  must  equate  the  gradient  of 
o2  to  zero.  This  procedure  determines  two  equations  for  the  mini¬ 
mum  point: 

I 

< ( R-r ) • i > - <n • ( R-r)  n*i>  = 

<(R-r)*j '  -  <n-(R-r)  n-j>  =  .  ( A- 6 ) 

Expressing  all  vectors  in  terms  of  their  cartesian  components , 
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n  "  1  Kind  +  j  c ok 4 
a  Y  a  d 


r 

~a 


/V 

i 


X 

a 


A 

+  j  y 


a 


K 


X  +  j  Y 
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we  obtain  the  following  pair  of  coupled  linear  equations: 


r  <eosz<J>>  -  <z  in^cos <f>>\ 

n. 

I  <cos<p(xcoc<p-  zin<p)>  ] 

i  -  <5in<j>cos <p>  <sin2<^>  ; 

iV  1 

{-  <cln4>(xcos0-ysin^)>  t 

The  condit  ion  for*  the  existence  of  solutions  of  £q,  A-7  Is 
that  the  determinant  of  the  coefficient  matrix  docs  not  vanish. 
Note  that  the  coefficient  matrix  Ls  entirely  independent  of  the 
array  locations;  the  array  locations  enter  only  in  the  inhomoge- 
neous  term.  v  me  a?  is  of  vlmpl*-*  algebraic  manipulations,  re  can 
express  the  d*-*t  erminan*  in  the  form 

DKT  =  -  cos?4i-<-cou?<t>>  )?^>  +  <^(sin?({--<.  in2<i>  )2  i  .  (A-8) 

The  determinant  Is  the  sum  of  the  mean-square  deviation  of  the 
cosine  of  twice  the  azimuth  and  the  sine  of  twice  the  azimuth. 

oth  of  these  .terms  can  vanis  only  if  all  of  the  azimuths  are 
equal.  ,\e  conclude  therefore  t..at  the  MD  prediction  fails  to 
exist  only  if  all  the  bearing  lines  are  parallel. 

If  the  determinant  does  not  vanish,  the  solution  of  Mq .  A-7 
for  the  coordinates  of  the  MD  prediction  is 
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4 


X 


<  u  I  n  *  ip><cosip(  x  ■  oo  <t>-yoin<f>)>-<rim}>coi,4>>*  a  •  in^(  xcooif  -  _  3in4>)> 
<o  in2  <}>><cos2<})>  -  [<cimf>cos<fi>  ]2 


Y 


<3in!|)GOS  (ftxcootj--  (  xco.'  c*  —  y s i n )  >-< 2  4>><  3  1  n<f> (  xcoo<{  -  sln4> )  > 
<o  in2  ({>><cos24>>  -  [  <3 in<t>CG3  <p>  ] 2 


(A-9) 


Having  obtained  valuer  of  X  and  Y  “rom  this  equation,  we  can  find 
the  rmo  deviation  of  the  prediction  from  the  bearing-  lines  directly 
with  Eq.  .  A- I  and  A -lj  . 
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APPENDIX  B 


Generation  of  Synthetic  Data  for  Testing  of 
the  Generalized  Tri an gu 1  at i on  Algorithms 

We  assume  a  four-array  configuration  of  azimuth-determining 
array.-  as  shown  in  Fig.  b-1.  Mach  individual  array  is  the  same  as 
In  Fig.  2.  We  postulate  a  group  of  19  targets  at  locations  shown 
on  the  map  of  Fig.  8-2.  The  range  and  azimuth  of  each  target  rel¬ 
ative  to  the  center  of  the  configuration  (located  at  the  central 
microphone  of  array  4)  is  given  in  Table  B-I.  The  target  field 
includes  three  closely  spaced  clusters  along  lines  making  differ- 
nt  angles  with  the  baseline  of  the  configuration. 

We  postulate  a  series  of  one  or  more  shots  originating  at  each 
target  at  times  given  in  Table  B-I .  Some  targets  are  responsible 
for  shots  spaced  less  than  thirty  seconds  apart.  In  addition,  the 
shots  emanating  from  each  cluster  of  targets  all  occur  in  the  same 
two-  or  three-second  interval.  The  shot  times  are  referred  to  the 
time  of  the  first  shot,  which  occurs  at  target  17*  We  number  the 
shots  from  one  to  th irty-se ven  in  order  of  their  times  of  initia¬ 
tion. 

From  the  geometry  of  the  target  and  array  locations,  we  cal¬ 
culate  the  true  azimuth  of  each  target  relative  to  the  center  of 
each  array  (Table  b-II ..  These  art  the  azimuths  which  would  be 
determined  by  a  perffcet  azimuth  determination  algorithm. 

.he  next  step  is  polluting  the  azimuthal  data  to  simulate  the 
results  of  an  irnperf  c*  azimuth  determination  algorithm  (which  is 
subject  to  azimuthal  Quantization  and  other  errors).  ,.’e  altered 
each  of  the  148  shot  azimuths  (four-  azimuths  for  each  of  37  shots) 
accoraing  to  the  following  scheme: 

(i)  To  each  azimuth,  add  an  increment  of  either  -0.75°,  -0.25°, 

+  0.25°,  or  +0.75°  with  a  one-fourth  ;  rebability  for  each  in¬ 
crement  . 
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EQU I L  ATEP AL  APR  A v  CONFIGURATION 


TABLE  B  - 1  I .  True  target  azimuths  at  each  array  (to  nearest  5'). 


Target 

Azimuth  at 
Array  1 

_ I _ 

Azimuth  at 
Array  2 

_ _ _ . 

Azimuth  at 
Array  3 

Azimuth  at 
Array  4 

1 

96° 

2C  ' 

36° 

30' 

93° 

51  ' 

42° 

10’ 

- 

39° 

0r  ' 

22° 

'  0’ 

33° 

11  ' 

•iO  0 

10' 

3 

72° 

' 

1  8° 

95' 

-  9° 

30 ' 

65° 

' 

357° 

31  ' 

350° 

3f  ' 

5° 

51  ’ 

356° 

' 

5 

0° 

0C  ’ 

352° 

20' 

7° 

90 ' 

O 

J 

00  ' 

t 

TO 

4_ 

lr  ' 

3  ‘ji\° 

00' 

9° 

25' 

1° 

50’ 

7 

14° 

9‘  ' 

358° 

55' 

11° 

15' 

9° 

00  • 

8 

7° 

? 

3  57° 

9 r  ’ 

13° 

' 

ltO 

91  ’ 

0° 

.  ' 

J  -J 

350° 

35’ 

19° 

50' 

8° 

’ 

10 

3J*2° 

or  t 

5  39° 

3r  ’ 

353° 

r  ' 

3')5° 

0  0  ' 

1  ) 

■A  i 

3^3° 

t 

390° 

’ 

353° 

-)  r  » 

395° 

1  o 

i.  c 

39  3° 

4ir  f 

3AO0 

31  ' 

35  3° 

30' 

39  5° 

51  ’ 

1 3 

■i  n  n  ° 

i  • 

3  A 1  0 

' 

353° 

9  0  ' 

39o° 

15’ 

m 

3li  ii 0 

3  ' 

•  9 1 0 

' 

353° 

50' 

396° 

90  ' 

t  c-i 
i  j 

302° 

? 

i  1 

ni° 

1  1  ’ 

32  >° 

30  ' 

319° 

50’ 

1  0 

30  i° 

V  ' 

313° 

11  ' 

'i  o  n 

<  t 

5' 

315° 

00’ 

■t  *7 

2  8  7  ° 

1 

2  9° 

2  5 » 

js25° 

30 ' 

306° 

00' 

i  8 

288° 

’ 

30  3° 

31  ' 

-i.  9° 

' 

309° 

50' 

l  cj 

288° 

' 

303° 

’ 

o  i  o 

3  a.  ji 

0’ 

309° 

0C  ' 

(ii)  Round  off  each  azimuth  to  the  nearest  half  degree. 

Tli.e  distribution  of  error.'  resulting  from  this  procedure  is  shown 
in  Fig.  B-3.  The  final  shot  azimuth  sets  for  all  37  shots  are 
presented  in  Table  This  is  the  data  we  use  to  test  the 

generalized  trianrul ation  methods  in  Sec.  3- 


NUMBER  OF  ERRORS 


TABLE.  B-III 


Shot  azimuth  sets  with  simulated  errors 


Azimuth  at 
Array  2 


30* .  50° 
30 A .  •''0° 
3n2 . 50° 
37.00° 
22.50° 
58.50° 
59.00° 
36.00° 
22.50° 
352.00° 
355.50° 
363.60° 
350.50° 
367 .  r'0° 
360.00° 
339.00° 
339-50° 
3*1.00° 
3*0.50° 
3*3*00° 
37. oo° 
22.50° 

313.50° 

59.00° 
312.50° 
22.50° 
37.00° 
5  9 .  00° 
22.50° 
37.00° 
37. 00° 
22.50° 


313. 

311.50* 


Azimuth  at 
Array  3 


326.00° 
32* .00° 
323. 00° 

*3-50° 
33.50° 
6*  .00° 
6* .00° 
** . 50° 
33.00° 
8.00° 
11.50° 
9.  0° 
5.50° 
13.50° 
1  * . 50° 
353.50° 
353.  0° 
353.00° 

3*3.  0° 

35* .00° 
** .  60° 
33.00° 
322.50° 
6*.  0° 
329.00° 
33-50° 
*3-50° 
r  5.00° 
33-50° 
*3-50° 
**.  0° 
33-50° 
33-  0° 
322.50° 
330.00° 
321.50° 
329.00° 


Azimuth  at 
Array  4 


305.50° 
305.60° 
303.50° 
*2.50° 
29.50° 
6*. 50° 
65.50° 
*2.00° 
30.50° 
.  50° 
3-60° 
1.50° 
358.50° 
5.50° 
7.50° 
3*5.50° 
3*5.50° 
3*5-50° 
3*6.50° 
3*7.00° 
*2.50° 
30.50° 
31*. 50° 
65.50° 
31*. 50° 
30.50° 
*1.50° 
65.50° 
29.50° 
*2.60° 
*2. r0° 
29-50° 
30.  .0° 
31*. 50° 
31* . 50° 
315.50° 
31* .50° 
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This  report  considers  the  prol  lem  of  predicting  target  coordinates, 
given  a  temporal  and  spatial  distribution  of  shots  from  a  multi* 
pllcity  of  targets.  ur  method  is  based  on  analysing  multifold 
polarity -coincidence  correlations  of  acoustic  signals  received  at 
a  configuration  of  asimuth- Determining  arrays. 

We  review  the  process  of  generating  multifold  correlation  functions, 
with  emphasis  on  those  properties  which  are  of  concern  in  azimuth 
determination.  We  describe  an  azimuth  determination  algorithm 
designed  to  make  most  effective  use  of  the  information  contained  in 
the  correlation  function. 


Our  discussion 
of  traditional 
tri angulation  t 
of  azimuths . 
demanding  in  co 
method  based  on 
This  new  method 
data,  and  found 
method . 


of  generalized  triangulation  points  out  the  weaknesses 
methods  ,  and  outlines  a  probabilistic  approach  to 
hat  would  yield  an  optimal  prediction  for  a  given  set 
ecause  this  probabilistic  approach  would  be  extremely 
mputational  requirements,  we  present  an  alternative 
similar  concepts  but  yielding  a  closed-form  solution, 
of  generalized  triangulation  ir  tested  with  synthetic 
consistently  more  accurate  than  the  traditional 
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